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) ^ ' We study time evolution of the Q ball in thermal logarithmic potential using lattice simulations. 

^— i , As the temperature decreases due to the cosmic expansion, the thermal logarithmic term in the 

r— ^ ■ potential is eventually overcome by a mass term, and we confirm that the Q ball transforms from 

*y\ the thick-wall type to the thin-wall type for a positive coefficient of radiative corrections to the 

mass term, as recently suggested. Moreover, we find that the Q ball finally "melts down" when the 
Q-ball solution disappears. We also discuss the effects of this phenomenon on the detectability of 
i-G ' gravitational waves from the Q-ball formation. 
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I. INTRODUCTION 



CN 

J> ■ A Q ball [1| is a nontopological soliton, which consists of scalar fields that carry global U(l) charge Q. Its existence 

I/"") | and stability are guaranteed by finite Q. Q balls are interesting objects in cosmology because they are often generated 

C*~) . in the Affleck-Dine (AD) mechanism for baryogenesis and can play an important role for baryogenesis and be a good 
dark matter candidate 043 ■ Recently, it was claimed that gravitational waves (GWs) are generated at the Q-ball 
formation |{|, which may be detected by the next generation gravitational wave detectors such as DECIGO [1(| and 

C~~- ■ BBO [llj. However, the detailed study of the subsequent evolution and the decay of Q balls revealed it to be difficult 

( ^ ) even by those next generation gravitational wave detectors [12( . 

The properties of the AD mechanism and the Q ball depend on the supersymmetry (SUSY) breaking mechanism 
@j 3 j because the effective potential of the relevant scalar field (AD field) quite differs for the different mediation 
mechanism. Consequently, there are various types of Q balls, such as gauge-mediation type [2|, |4[ , gravity-mediation 
type @, Q , new type @ , delayed type @, Q , and so on. Among them, the thermal log type Q ball @ , whose effective 
potential is dominated by the thermal logarithmic term [l3[, possesses an interesting feature. As the Universe expands, 
the cosmic temperature decreases and so does the thermal logarithmic potential. Thus the properties of thermal log 
type Q ball will change with time. Moreover, zero-temperature potential eventually overcomes the thermal potential 
and then the type of Q ball changes [12( . If the zero-temperature potential itself allows a Q-ball solution, the type of 
Q ball changes to the corresponding type. 

It may then be naively expected that the Q balls would be destroyed if the zero-temperature potential alone does not 
allow a Q-ball solution. Recently, however, one of the present authors (S.K.) showed that even if the zero-temperature 
potential alone does not allow a Q-ball solution, the total potential (the thermal logarithmic term and a mass term 
with a positive radiative correction) does allow a Q-ball solution, which would result in the transformation from the 



thermal log type of the Q ball into the thin- wall type [14|. Since the scenario would be changed in this case, it is 



important to investigate whether the field configuration dynamically transforms from one type of the Q ball to the 
other. 

In this article we perform numerical simulations on the lattice to see the time evolution of the configuration of the 
AD field in the potential with a thermal logarithmic term and a mass term with a positive coefficient for radiative 
corrections, where the latter term alone does not allow a Q-ball solution. We confirm that the thermal log type Q 
ball transforms to the thin- wall type Q ball found in Ref. [Tj]. In addition, we find that the Q-ball configuration 
"melts down" when the cosmic temperature becomes too low to hold a Q-ball solution. Wc also find that it changes 
the thermal history of the Universe, but the detectability of the GWs from the Q-ball formation is not improved 
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compared to the estimate of our previous study of Ref. [12J . 

The paper is organized as follows. In the following section, we review the properties of the thermal log type Q ball 
and of the Q-ball solution found in Ref. [14[, which we will call the thermal thin- wall type. In Sec. IIII[ we show the 
numerical results of the time evolution of the field configuration and confirm the transformation of the Q-ball types. 
In Sec. IIV1 we see numerically the "melting down" of the Q ball as the temperature decreases further. In Sec. [V] 
we reconsider the fate of the Q ball and the AD field in this case. We also discuss the effect of this feature on the 
dctectability of the GWs from the Q-ball formation. Finally, Sec. IVII is devoted to the conclusion. 

II. PROPERTIES OF Q BALLS 

We are interested in the Q-ball properties in the potential with both a thermal logarithmic term and a mass term 
with a positive coefficient for one-loop radiative correction, 

(1) 

(2) 

(3) 

where $ is the complex AD field and T is the cosmic temperature. The upper term in Fthormai represents the 
thermal mass from the thermal plasma and the lower one represents the two-loop finite temperature effects coming 
from the running of the gauge coupling g(T) which depends on the AD field value [l^. 1 Note that even before the 
reheating from the inflaton decay has not completed, there exists thermal plasma from the partial inflaton decay as 
a subdominant component of the Universe. Ugrav denotes a soft mass term due to gravity-mediated SUSY breaking, 
where m^ ~ O(TcV). The second term in the bracket is the one- loop radiative correction, and A is the renormalization 
scale. Here we assume K > (K ~ 0.01 — 0.1) so that Vg rav alone does not allow a Q-ball solution. 

At larger temperature when the field starts the oscillation, the potential is dominated by the thermal logarithmic 
term Vthermab and the thermal log type Q balls form @, Uj]. The properties of this type Q ball are similar to those 

MT)~TQ l /\ W (T) ~ V27TTQ- 1 / 4 , E(T) ~ ^TQ 3 / 4 , R(T) ~ ^, (4) 

where Q is the charge stored in a Q ball, <pa = v / 2| < I , o| is the AD field value at the center of Q ball, u is the angular 
velocity of the AD field, E is the energy stored in a Q ball, and R is its radius. Since the charge Q is the conserved 
quantity, whose value is determined at the Q-ball formation, the parameters of Q balls change as the temperature 
decreases according to Eq. (UJ). Notice that the configuration of the AD field will follow the Q-ball solution, since the 
time scale of the Q-ball reconfiguration is much shorter than the cosmic time: T _1 <C H -1 . 

As the temperature decreases further, V gTav will eventually dominate the potential at <f>Q. In our previous study in 
Ref. [T2J, we assumed that Q balls are destroyed and turn into almost homogeneous AD field quickly at this moment, 
because the potential ([3]) alone does not allow a Q-ball solution. Recently, however, one of the present authors pointed 
out that a Q-ball solution does exist even in this situation [14| . Although the soft mass term overcomes the thermal 
logarithmic term at large field values, the latter will dominate the potential at smaller field values. As a result, in the 
light of charge conservation, a thin-wall type Q-ball solution exists. 

Let us investigate the condition of the existence of this Q-ball solution. A Q-ball solution exists when V/(f> 2 has a 
global minimum at </> = <^ m in7^0 |l|. In this case, it is satisfied when 
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has a nonzero real solution. In fact, a Q-ball solution exists for T > m^K 1 / 2 , even after the potential is dominated 
by the soft mass term at the center of the Q ball. Wc shall call it the thermal thin- wall type Q ball. 



We neglect the numerical factors coming from the gauge coupling constants. 



Now we study the properties of this type of Q ball further. A thin-wall type Q-ball solution satisfies [lj 

^L R 3 ~ Q (6) 

B~ Q f^. (T, 

For given Q, (f>o is determined so that the energy of the Q ball should be the minimum. Thus the properties of the Q 
ball arc written as 

M?) ~ c(T/n H K^ 2 ) T ej~a(T)m*, E ~ a(T)m^Q, R~ (^f£) , (8) 

m<f,K L i z \ cal 4 J 

where c(T /m^K 1 ' 2 ) and a(T) are slowly increasing functions of T and they are of order of unity at the temperature 
in which we are interested. For example, c(10) ~ 2.5,c(10 2 ) ~ 3.4,c(10 3 ) ~ 4.1,c(10 4 ) ~ 4.6,c(10 5 ) ~ 5.1,c(10 6 ) ~ 
5.5,c(10 7 ) ~ 6.0 and so on. a(T) is expressed as 

a2 = 1+K ( log (2^) + h log {^k) ) ' (9) 

and its temperature dependence is stronger than that of c. 

It is true that such a Q-ball solution exists but not clear that the field configuration follows from the thermal log 
type Q ball to the thermal thin-wall one. Moreover, it is nontrivial what happens when the Q ball solution vanishes 
at the temperature T ~ ra^K 1 / 2 . In order to tell how the configuration of the AD field evolves, we perform numerical 
studies on the lattices for each case in the following two sections. 

III. TIME EVOLUTION AND TYPE TRANSFORMATION OF THE Q BALL 

In this section we investigate the evolution of the Q ball in the potential of the thermal logarithmic term and the 
soft mass term with positive radiative corrections. Since we arc primarily interested in the transformation of a Q 
ball, here we limit ourselves to a single Q ball assuming the spherical symmetry of the field configuration, and solve 
the one-dimensional partial differential equations in the radial direction by using the staggered leapfrog method with 
second order accuracy both in time and in space. 

In order to obtain an initial configuration of the field on the lattices, we first solve the ordinary differential equation, 

ar z r dr \ a<p J 

with boundary conditions, d<p/dr{r = 0) = and <p(r = 00) = by using the fourth order Runge-Kutta method. 
Instead of Eq. @ , here we use the following thermal potential: 

2T 2 



F ther mal = T 4 l0g 1 + ^3 , (11) 



in order to interpolate the thermal mass term at smaller field and the thermal logarithmic term at larger field values. 
This form of potential includes the two limit in Eq. ([2]) and connects them smoothly. Rescaling variables with respect 
to A [in Eq. ((3}], we use the following dimcnsionless variables: 

4>=j, r = rk, E = T i=M ' ^ = T (12) 

Since we are seeking for the initial configuration when the thermal logarithmic term dominates the potential where 
the Q ball is the thermal log type, we set T*/A = 1.0 x 10 -2 as the initial temperature and m^/A = 1.0 x 10 -4 . 
Other parameters are set to be K = 0.1 and Mq/A = 24.3, where Mq is the reduced Planck mass. Configurations 
which we find are summarized in Table HI They coincide with the evaluation in Eq. Q within a numerical factor of 
order of unity. 

Now we can investigate the time evolution of the Q ball by using the second order leapfrog method. Hereafter 
we assume the radiation dominated universe. In order to guarantee the regularity at the origin, we use the variable 



TABLE I: Properties of the thermal log type Q balls. 



Charge (Q) Angular velocity (o>) 


Energy (E) 


Field value (<f>o) Radius (R) 


2.5 x 10 8 6.0 x KT 4 


3.9 x 10 5 


2.9 3.2 x 10 3 


1.0 x 10 9 4.4 x 10~ 4 


1.1 x 10 6 


4.2 4.5 x 10 3 


3.0 x 10 9 3.5 x 10~ 4 


2.7 x 10 6 


5.5 6.0 x 10 3 



X = r/fcj). r\ (= 2-y i*i '/ 'a*) is the conformal time rescaled with respect to A, where t* = [90 / (n 2 g*)] 1 ^ 2 AMq / (2T 2 ) 
and g* ~ 200. a* is the scale factor at the temperature T = T* and set to be 1. In the radiation dominated universe, 
T oc rj (x a -1 . We use the time step as dr\ = 0.5. The grid spacing is df = 1.8 x 10 5 /2 16 ~ 2.7 and the number of grid 
is 2 16 . We decompose the real and imaginary parts of the scalar field as xe lut = xr + ix.1- I n terms of r\ and r, the 
evolution equation of Xi {i = ^ -0 becomes 



d Xi d 2 Xz , a 2 r]f dV tot 



drj 2 



df 2 



A 4 d Xl 



0. 



The boundary conditions at the origin (f = 0) are set to be 



Xi(f = o) = 0, 



dXi 

drj 



(f = 0) = (i = R,I), 



(13) 



(14) 



and we assign the free boundary conditions at the other end of the grid. We have checked the charge conservation 
both in an expanding and nonexpanding universe and the energy conservation in a noncxpanding universe and have 
confirmed that all hold with an accuracy of 10~ 7 throughout the calculation. Thus, we conclude that our numerical 
simulation is accurate enough for our study. 

The time evolution of the field configuration is shown in Fig. [T] for Q ~ 1.0 x 10 9 . The axes are rescaled with 
respect to the scale factor a so that the rescaled radius is almost constant for the thick-wall type. We can see that 
the configuration of the Q ball changes from the thick-wall to the thin-wall types. This coincides with the feature of 
the transformation of the Q-ball solution found in Ref. [1J| . 
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FIG. 1: Configurations of the AD field for Q ~ 1.0 x 10 9 at the time of T = T„, T»/2, T»/4, T,/8, T./16, and T./32 from the 
top to the bottom, respectively. Q-ball configuration changes from the thick-wall to the thin-wall types. 



Figure [2] shows the temperature dependence or the time evolution of the rotation, gradient, and potential energies 
of the field configuration for Q ~ 1.0 x 10 9 . The gradient energy decreases with temperature more quickly than other 
contributions, and both the rotation and the potential energy decrease in a similar manner. These features also imply 
that the thermal log type Q ball becomes flattened and behaves like the thin-wall type one. 

So far we see the transformation of the Q-ball configuration: the thick-wall to the thin-wall type. We further check 
the properties of the Q ball before and after the transformation to confirm that they are those of the thermal log 
type and the thermal thin- wall type, respectively. The temperature dependence of the field value at the Q-ball center, 
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FIG. 2: Temperature dependence of the gradient, rotation, potential, and total energies of the AD field. 



the Q-ball radius, 2 and the field angular velocity with various charges of the Q ball are shown in Fig. [3] Here (blue) 
crosses, (red) x's, and (green) stars represent numerical results for Q ~ 2.5 x 10 8 , 1.0 x 10 9 , and 3.0 x 10 9 , respectively. 
Corresponding lines are the analytic estimates ^ and (JSJ up to numerical coefficients for each case (short dashed 
lines (blue): Q ~ 2.5 x 10 8 , straight lines (red): Q ~ 1.0 X 10 9 , long dashed lines (green): Q ~ 3.0 x 10 9 ): At high 
temperature, <po °< TQ X / A ,R ex T~ 1 Q 1 ^ 4 ,oj ex TQ -1 / 4 for the thermal log type Q ball, while, at low temperature, 
</>o oc T 2 ,R ex T -4 / 3 ^ 1 / 3 ,^ = const, for the thermal thin-wall type Q ball. We show the analytical estimates of 



the field value and the angular velocity in Eq. <j8j> with dotted lines (pueplc) in Fig. 3(a) and 3(c) since they are 
independent of charge Q. We can see that the analytical estimates ((4]) and ([SJ are well reproduced by the lattice 
simulations. One exception is the angular velocity at low temperature. This could be understood by the factor a(T) 
in Eq. ([8]), which decreases as T gets lower. We can therefore conclude that the Q ball really transforms from the 
thermal log type to the thermal thin-wall type in the potential considered here. 

One interesting feature is the increase in the radius of the thermal thin-wall Q ball as the temperature goes down. 
Since the growth rate of the radius is larger than that of the cosmic expansion, Q balls may merge eventually. 



IV. MELTING OF THE Q BALL 



In this section we investigate the evolution of the field configuration around the time when the potential can no 
longer hold a Q-ball solution, T ~ ra^K 1 / 2 by using the numerical calculations based on the second order leapfrog 
method. To this end, here we set T*/A = 1.0 x 10~ 2 , m^/A = 1.0 x 10~ 2 , M G /K = 24.3, and K = 0.01 to seek the 
thin-wall type Q ball as an initial configuration of the field on the lattices by using the fourth order Runge-Kutta 
method, and find the configurations in Table HT1 They coincide with the evaluation in Eq. (J8j) within a numerical factor 
of order of unity. We then solve the equations of motion Eq. (fT3| to see the time evolution of the field configuration 
on the lattices. Here we use the grid spacing df = 4.1 x 10 6 /2 19 ~ 7.8 and the number of the grid 2 19 . We have 
confirmed the charge conservation is held with an accuracy of 10 . 



TABLE II: Properties of the thermal thin-wall type Q balls. 
Charge (Q) Angular velocity (w) Field value (<j!>o) Radius (R) 



4.0 x 10 9 

1.1 x 10 10 
3.0 x 10 10 



1.0 x 10" 2 
1.0 x 10 -2 
1.0 x 10~ 2 



0.56 
0.53 
0.52 



9.3 x 10" 
1.3 x 10 4 

1.8 x 10 4 



Figure |4] shows how the AD field configuration breaks down in the case Q ~ 1.1 x 10 10 . Here the vertical and 
horizontal axes are rescaled with respect to a 2 and a~ 4 / 3 , respectively, so that the rescaled radius and field value are 
almost constant for the thermal thin-wall type. Temperature dependence of the radius of the Q ball or AD field lump 



We define the radius as the length between the center of the Q ball and the point where <j> = </>o/2. 
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(b) Radius of Q ball 
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(c)Angular velocity of the AD field 

FIG. 3: Temperature dependence of Q ball properties. Crosses (blue), x's (red), and stars (green) represent the numerical 
results for Q ~ 2.5 x 10 8 , 1.0 x 10 9 , and 3.0 x 10 9 , respectively. Lines correspond to analytical estimates (UJ and J8| up to 
numerical coefficients. 

is shown in Fig. where the radius is defined as the length between the center of the Q ball and the point where 
<f> = </>o/10. Crosses (blue), x's (red) and stars (green) represent the numerical results for Q ~ 4.0 x 10 9 , 1.1 x 10 10 , and 
3.0 x 10 10 , respectively. Corresponding lines are the analytic estimate (j8]) at high temperature (R oc T~ 4 ' 3 oc a 4 ' 3 ) 
up to numerical coefficients for each case (short dashed line (blue): Q ~ 4.0 x 10 9 , straight line (red): Q ~ 1.1 x 10 10 , 
long dashed line (green): Q ~ 3.0 x 10 10 ). Dotted line (purple) shows the growth rate expanding at the speed of light 
(R ~ i*(T/T*) -2 ) up to the numerical coefficient. We can see that the field configuration follows the thin- wall type 
of the Q-ball solution at high temperature. At the temperature T ~ m^K 1 / 2 , where T ~ 0.1T* for the parameters 
taken in Figs. 0] and [5J the Q-ball solution vanishes and the wall of the Q ball starts "melting down." At lower 
temperature, the AD field lump expands at the speed of light and the configuration becomes homogeneous. This 
behavior is consistent with the condition, Eq. ([5]). 

Therefore, we can conclude that the Q ball starts "melting down" when the potential can no longer allow a Q-ball 
solution and its remnant spreads out at the speed of light, leading to the homogeneous AD field configuration in a 
few Hubble time. Note that the Q ball may have already decayed before the "melting" time, since, at that time, the 
field value at the center of the Q ball becomes </>o ~ m^K 1 ' 2 , which is as small as the electroweak scale, O(TeV). 
At that time the decay channel to the fields coupled to the AD field opens and it would decay before or during 
homogenization. 



V. LATE TIME BEHAVIOR OF THE Q BALL 



In the previous sections we confirmed that Q balls in the potential with the thermal logarithmic term and the mass 
term that receives a positive radiative correction transform from the thick-wall type to the thin-wall type and then 
"melt down" into the homogeneous AD field. In this section, we consider the late time behavior of the time-dependent 
Q ball and its decay time in detail. 

As opposed to the thermal log type Q balls whose energy density decreases as T 4 , the thermal thin- wall type Q 
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FIG. 4: Configurations of the AD field for Q - 1.1 x 10 10 at the time of T = T*,T»/2, T»/4,T»/8, T./16, and T./32. 
configuration starts "melting down" and expands. 
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FIG. 5: Temperature dependence of Q-ball radius. Crosses (blue), x's (red) and stars (green) represent the numerical results 
for Q ~ 4.0 x 10 9 , 1.1 x 10 10 and 3.0 x 10 10 , respectively. Corresponding lines are the analytical estimates (O up to numerical 
coefficients. Dotted line (purple) represents the line R ~ £»(T/T«)~ 2 
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balls and the homogeneous AD field behave like matter. Therefore, they would eventually dominate the Universe, 
which would alter the estimation of the detectability of the GWs from the Q-ball formation in Ref. [12| . In order to 



investigate whether they really affect the cosmic history, we need to estimate the lifetime of the Q ball. 

There are two cases of the cosmic history: The Q balls will decay or merge. In order to see which will take place, 
we need to estimate the Hubble parameter at the Q-ball decay, Hd cc , and the temperature at the Q-ball merger, T m . 
Note that the fate of the Q balls differs for the different SUSY breaking mechanisms as discussed in Ref. [12[ . 

Let us first consider the gravity-mediated SUSY breaking mechanism. The decay rate of the Q ball can be estimated 
as follows. When the field value inside the Q ball is larger than the electroweak scale E w , the decay inside the Q ball 
is kinematically forbidden because the would-be decay products acquire too a heavy mass, g4>o, for the AD field to 
decay into, and hence the Q ball can decay only from its surface. Here g is the Yukawa/gauge coupling constant of 



order of unity. This sets the upper bound on the decay rate of the Q ball [1 

dQ 



1 

■1 dec, surf = TT 



at 



u 3 R 2 , N 

* 48^- ^ 



In fact, it is known to be almost saturated for u>/g4>o < 1 [Tf|. Thus the Q ball decays when the Hubble parameter 
becomes smaller than this decay rate if the field value in the Q ball is large enough. On the other hand, when the 
field value inside the Q ball gets smaller than the electroweak scale, T < (m^E^) 1 ^ 2 , or H < m^E^/Mc, the decay 
channel opens, g4>a/oJ < 1, and the Q ball can decay from the interior. At that time, the decay rate is estimated as 
in the case of homogeneous oscillating fields, T ~ h 2 m/(8Tr), which is much larger than the Hubble parameter. Thus 
the Q ball decays immediately. Therefore, the Hubble parameter at Q-ball decay is expressed as 




^ =max ^- f (487r)3/7c2/7gl/7 Mr w -McT/- («) 

In this estimate, we assume that the reheating temperature is high enough that the Q-ball decay takes place during 
the radiation dominated era. 

In the above argument, we have not taken into account the fact that the thermal thin- wall type Q balls would 
merge together and turn into an almost homogeneous AD field, 4 as mentioned in Sec. IIIII This would take place if 
H m>er > i?doc, S urf, where 

m ' gr " M G ' m ' gr - T R M% ■ [ l} 

Here 4> osc is the AD field value at the onset of Q-ball formation and Tr is the reheating temperature from the inflaton 
decay, and the distance between Q balls is estimated as 

L = L f x^ = L f xf^V /3 x— . (18) 

at \H R J T' V ' 

where Hr is the Hubble parameter at the reheating. The distance between the Q balls, the Hubble parameter and 
the Q-ball charge at the formation of the thermal log type Q ball during inflaton oscillation dominated universe are 
given respectively by [12|, [l4[ 



L { ~Hf\ Hf~lO- 2 -^, Q ^2xl0-^. (19) 

'Pose 1V1 G 1 R. 

The energy density of the almost homogeneous AD field is a little smaller than the total energy density of the Universe 
at that time: /3ad(— m 00cq — ^m) ~ Ptot{— (.9* 7I ' 2 /90)T 1 4 1 ), where cq ~ T^/m^. Thus it will soon dominate the 
energy density of the Universe. Note that the thermal thin-wall type Q balls will merge before they would dominate 
the energy density of the universe. After that, the AD field value becomes as small as the electroweak scale when 
H ~ m f pE w /Mc- At that time, the AD field decay is kinematically allowed and takes place immediately since the 
Hubble parameter is much smaller than the decay rate. This is the same as the scenario that the transformation of 
the Q ball into the thermal thin- wall type had not been taken account of, as in Ref. |12| . The only difference in the 
presence of the thermal thin- wall type Q balls is that they may decay earlier from their surface before they dominate 
the energy density of the Universe. 

Next we consider the Q ball in the gauge-mediated SUSY breaking mechanism. The potential is then given by 

Vtot = ^gaugc + V / grav2 + Khermalj (20) 

^ g au g c=M 4 fl0g^") , (21) 
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(22) 
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Originally, this bound is for the decay into fermions, which is forbidden in the interior of the Q ball due to the Pauli principle. However, 
we can use this bound for the decay into scalars since this bound is determined by the surface region of Q ball but not by the type of 
the decay product. 

The phases and the field values of the AD field in the Q balls differ with each other and hence Q balls may not merge easily. However, 
the merger may be completed in a few Hubble time, and hence we may assume that the merger takes place instantaneously at T migr . 



Here M$ is the messenger mass, Mp ~ (?7i^Ms) 1//2 (> O(lTeV)), and m 3 / 2 (< O(10 GeV)) is the gravitino mass, 
which is smaller than the electroweak scale. We are interested in the K > case where the thermal log type Q 
ball transforms into the thermal thin- wall type due to Vg raV 2. Once this transformation occurs, the thermal thin- wall 
Q balls will end up with one of the following three destinies. The first one is that they will merge together at the 
temperature 



Km 3/2 (j)t 
T R M% 



Jm,ga— t-, „, r3 , [to) 



and turn into a homogeneous AD field rotating in Vg raV 2. Since the field amplitude decreases as <j> oc a~ 3 / 2 oc T 3 / 2 
while cj) q oc T 2 where Thermal = U graV 2 at = oq , (/) never reaches to Vthcrmai- However, since Vtot is eventually 
overcome by Vg augc at (j> as the AD field value drops, the field <f> get s to Vg aug e to feel instabilities and the "delayed 
type" Q balls form. They will decay when H = r dccdo i ay , where [12j 5 

7T 2 m 3/2 
r d ec,dela y * ^^ M 4 • (24) 

This is the same as the scenario which we did not consider the transformation of the Q ball into the thermal thin-wall 
type in Ref. [l||. 

The second is that the thermal thin- wall type Q ball decays from its surface. Using Eq. (|!5p . we obtain the decay 
rate as 

(ir 2 9 A 2/J aK 2 / 7 miy 2 7 

J- dec, thermal — I 777^ I '. '. '. TT7 • V *-v 

V 90 ) (48^) 3 / 7 c 2 / 7 Q 1 /7 M 4/7 

The third is that the thermal thin- wall type Q ball ceases reconfiguration according to Eq.(J5} at T ~ Mp when 
Vgauge starts to overcome Vthermai- Henceforth, the thin- wall type Q ball has fixed configuration and decays when 
H = r dcc>fixod , where 

a 7/3 X 2/3 m ll/3 

Tdec, fixed = ; ; ^77 ■ (26) 

487TC 2 /3gl/3M^ /3 

Therefore, the transformation into the thermal thin-wall type could affect the cosmic scenario in such a way that the 
decay of the Q ball may take place later because the charge of the thermal and/or fixed thin- wall type Q ball will be 
larger than that of the delayed type. 

Finally, we comment on the influence of the thermal thin-wall type Q balls on the detectability of the GWs from 
the Q-ball formation. The present density parameter and the frequency of the GWs from the Q-ball formation can 
be generally written as [12[ 

*-*(!)■ (28) 

where the subscript "0" denotes that the quantity is evaluated at present. These values would change if there is 
an additional Q-ball dominated era. In our previous study [12|, we found that when Tr ~ 10 10 GeV, 4> osc — Mg, 
Mp ~ 10 4 GeV, and m 3 / 2 — 10 GeV in the gauge-mediated SUSY breaking mechanism, GWs from the Q-ball 
formation can be detected by DECIGO or BBO, but otherwise cannot be detected since their amplitude is too low 
and/or their frequency is too large. With these parameters, Q-ball merger takes place so that the cosmic history 
does not change. Thus, the detectability of the GWs will be the same. What about the cases with other parameters? 
In our previous study in Ref. [12| , we have found that the difficulties of the GW detection come from not only the 



5 If the mass of the decay product is larger than lo, the Q-ball decay is kinematically forbidden. In the case of the B ball that carries 
baryonic charge, u must be larger than the proton mass for its decay, which is hardly satisfied in the gauge-mediated SUSY breaking 
mechanism. However, in the case of L ball that carries leptonic charge but not baryonic charge, Q balls can decay if u) is larger than the 
lightest neutrino mass. In addition, in the case where the Q balls consist of a fiat direction with B — L = 0, there is no such constraint. 
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smallncss of the present amplitude of the GWs (^q W ^ 10~ 16 ) but also the largeness of the typical frequency at 
present (/o > 10 3 Hz). In the case of the gravity-mediated SUSY breaking mechanism, we see that the Q ball can 
decay earlier than the estimate in Ref. [12 |. which leads to a larger frequency of GWs. Thus the detectability of the 
GWs gets worse. In the case of the gauge-mediated SUSY breaking mechanism, we find that the decay rate of the Q 
ball may be much smaller than the estimate in Ref. [12J, which, in turn, leads to a smaller amplitude of GWs. Thus 
the detectability of the GWs gets worse also in this case. Therefore, taking into account the transition of the Q balls 
from the thermal log type to the thermal thin-wall type worsens rather than improves the detectability of the GWs. 

VI. CONCLUSION 

In this article, we have investigated the time evolution of the Q balls with both thermal logarithmic potential and 
the mass term with the positive radiative corrections. By means of one-dimensional lattice simulations for the radial 
direction, assuming a spherically symmetric profile, we have confirmed that the thermal log type Q balls transform 
to the thermal thin- wall type when the mass term overwhelms the thermal logarithmic potential at the center of the 
Q ball. In particular, we have found that the field value at the Q-ball center and the angular velocity of the AD field 
becomes the same value regardless of the charge stored in a Q ball. Furthermore, we have confirmed that the Q-ball 
configuration "melts down" when the cosmic temperature becomes too small to hold any Q-ball solution. 

Using this result, we have investigated the cosmic history in detail. Once the Q ball transforms into the thermal 
thin-wall type, the growth rate of the Q-ball radius becomes larger than that of the cosmic expansion. As a result, Q 
balls can merge before they decay and turn into almost homogeneous AD fields. Thus there are two scenarios in the 
fate of the Q balls: Q balls will merge together, or decay. If Q-ball merger takes place, the almost homogeneous AD 
field will dominate the energy density of the Universe soon after the merger. In the case that the Q balls decay before 
the merger, we have found that, for some parameter regions, the decay of the Q balls can be earlier in the gravity- 
mediated SUSY breaking mechanism and later in the gauge-mediated SUSY breaking mechanism, compared to the 
case without considering the transformation into the thermal thin- wall type, as in our previous study in Ref. [12J . 
However, the detectability of the GWs from the Q balls is not improved. 

Although we have limited ourselves to the spherical symmetric system and have studied the evolution of a single 
Q ball, which we think will capture most of the relevant features, it would be interesting to investigate the merger 
of multiple Q balls that may induce cosmologically interesting phenomenon such as another GW emissions. Since 
three-dimensional lattice simulations are required in order to see this in detail, it will be left as a future study. 
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